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Abstract 

It is shown that there exist symmetry constraints for non-leptonic weak am- 
plitudes which emerge when the 1/iVe-expansion restricted to the leading and 
next-to-leading approximations only is systematically combined with lim- 
ited to the lowest non-trivial order. We discuss these constraints for the couplings 
g 8 and g 27 of AS* = 1 transitions and the B^-parameter of K° — K° mixing. 
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1. Chiral perturbation theory (%PT) provides a useful framework to study K— 
decays in the Standard Model. To lowest order in the chiral expansion: 0(p 2 ), there 
are only two coupling constants g 8 and g 27 which govern non-leptonic i^-decays. Once 
these couplings are fixed phenomenologically from K —>■ tttt decays, there follow a 
wealth of predictions for K — > -r-r-k and some other radiative -fT-decaysf]. Unfortunately, 
the study of chiral corrections to these lowest order predictions generally brings many 
new local terms in the chiral effective Lagrangian of 0(p 4 ) 0-0. The number of 
terms is too large to make a systematic phenomenological determination of the new 
couplings similar to what has been done in the purely strong interaction sector ||. The 
problem here is not the complexity of the calculations; it is simply that the experimental 
information we have on AS = 1 and AS = 2 transitions is too limited when compared to 
the large number of possible C(p 4 ) couplings. Except for a few remarkable predictions 
which have been made, one is obliged in most cases to resort to chiral power counting 
arguments and/or models in order to make numerical estimates. 

Ideally, one would like to develop well controlled approximation methods starting 
at the level of the Standard Model Lagrangian. In that sense the l/iV c -expansion in 
QCD |9|, [TOfl, where N c — > oo with a s x N c fixed, is a good candidate. Keeping only 
the leading contributions in the large-iV c limit for non-leptonic i^-decays is however 
a bad approximation because, in that limit, many of the four-quark operators of the 
effective Lagrangian which emerge after integrating out the fields of heavy particles in 
the presence of gluon interactions are suppressed. One has to go to the next-to-leading 
order in the 1 /iV c -expansion before the complete set of possible four-quark operators 
appears. 

The purpose of this note is to show that there exist symmetry constraints for non- 
leptonic weak amplitudes which emerge when the 1/A^-expansion, restricted to the 
leading and next-to-leading approximations only, is systematically combined with xPT 
at the lowest non-trivial order. Here we shall limit ourselves to spell out these con- 
straints for the couplings g 8 and g 27 and the B^-parameter of K° — K° mixing and to 
the discussion of their phenomenological implications. There are other interesting ap- 
plications of the same type for other processes; in particular for the decay — ► 7r°e + e~ 
which at this approximation can be calculated in terms of known physical parameters, 
and which of course has interesting implications for K^ L — > 7r°e + e - and the possibil- 
ity of observing direct CP-violation in this process. These other applications will be 
discussed elsewhere. 

2. In the conventional formulation of xPT the octet of low-lying pseudoscalar states 
(-7T, K, rj) are the Nambu-Goldstone bosons associated to the "broken" axial generators 
of chiral-S'[/(3). The Nambu-Goldstone fields are collected in a unitary 3x3 matrix 
U(x) with det U = 1, which under SU(3)z £g> SU(3)r transformations (Vl,Vr) trans- 
forms linearly: U — > VrUV^. In order to describe non-leptonic weak interactions it is 
useful to introduce the 3x3 flavour matrix vector field 

*For recent reviews see e.g. refs. || 
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-i^U(x)^D fl U(x) 



where denotes the covariant derivative in the presence of external SU(3)l and 
SU(3)r gauge field sources, and / the /^-coupling in the chiral limit (/ ~ 86 MeV). 
Under chiral-5'C/(3) transformations: £ M — ► VlC/jVl- In terms of £ M , and to lowest 
order in the chiral expansion, the operators with the same chiral transformation prop- 
erties as those of the effective four-quark Lagrangian can then be readily obtained. 
They are: 



c 8 (x) = YtcMw) 



(2) 



and 



C 27 (x)^-(C^ 21 (^) n + (C^(D 



111 



(3) 



which transform respectively like (8l, 1r) and (27 l, under SU(3)l <%> SU(3)r. To 
lowest order in the chiral expansion, the effective Lagrangian of the Standard Model 
which describes AS = 1 transitions between pseudoscalar states has then the following 
form: 



r AS=l 
'"eft 



G 

KdKTs [g 8 A + g 27 ^27] + h-c 



(4) 



with g 8 and g 27 coupling constants which are not fixed by chiral symmetry arguments. 
The phenomeno logical determination of these couplings from K — > ttti, to lowest order 
in the chiral expansion, givesQ [[TTf 



|g 8 + gg 



3271 



5.1 



3271 



0.29. 



(5) 



The decays K — > and K — > 777777 have also been analyzed in the presence of chiral 
0(p 4 ) corrections ||12j| . The fitted value for g g decreases then by 30% to g 8 ~ 3.6 while 
g 27 is only slightly modified. 

It is useful to go one step backwards in the theory and to analyze the combinatorics 
which in the Standard Model leads to the effective Lagrangian above. With one virtual 
Vy-field emitted and reabsorbed, and to lowest order in the chiral expansion, there are 
three possible chiral invariant configurations which give rise to the effective Lagrangian 
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^Notice that £27 generates both AI — 1/2 and AI = 3/2 transitions. 
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where Q L are the flavour matrices 
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which under chiral rotations transform like — > VlQ^Vl- The underlying func- 
tional integral over quark and gluon fields which gives rise to the effective couplings 
in (P) is represented diagrammatically in Fig. 1. The solid lines correspond to quark 
fields propagating in a gluon background (the dots in the figure,) which is subsequently 
integrated down to the scales where the chiral Lagrangian of the Goldstone modes 
becomes effective. The Q^-operators represent the emission and absorption of the 
virtual W^-field. The restriction to lowest order in the chiral expansion implies that at 
most two insertions are allowed. When further restricted to AS* = 1 transitions, 
the effective Lagrangian in coincides with the conventional one in ([|) with 

3 3 
g 8 = -(a + b) -b + c and g 27 = -(a + b). (8) 



3. Let us now examine the behaviour of the coupling constants a, b, and c from 
the point of view of the l/iV c -expansion. It appears that the configuration which leads 
to the a-type coupling in (|]) is 0(N%), while those leading to the b and c type 
couplings are non-leading 0(N C ). To leading order in the l/iV c -expansion the coupling 
a can be calculated because in this limit the four-quark operators factorize into current 
density operators and their chiral effective realization is known from low-energy strong 
interaction physics to 0(p A ). With the factor f 2 , which is 0(N C ), included in the 
definition of C^(x) in (|lj), the coupling constant a is of 0(1) in the l/iV c -expansion. 
The interesting observation is that the factorization result which emerges, a = 1, can 
only be modified by gluonic configurations which are at least next-to-next-to-leading 
order in the l/A^-expansion, as illustrated by the diagram in Fig. 2. Colour matrices 
are traceless, which implies that a minimum of two gluons exchanged from one fermion 
loop to the other are required to modify the factorization property, and this leads to a 
relative correction of 0(1 /N 2 ). We then conclude that: 

a=l + o(^). (9) 

The configuration which leads to the c-type coupling corresponds to the so called 
penguin-like diagrams in the effective four-quark Hamiltonian formulation. It is well 
known ]13[ that their contribution to the coupling c to leading order in the 1/N C - 
expansion which in this case is 0(1 /N c ) can also be calculated in terms of known 
phenomenological parameters, with the result^]: 

$A detailed discussion of this calculation can be found e.g. in ref. 0] 
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c = ReC 4 -16L 5 ReC 6 (> 2 



<ipip> 



0.3 ±0.2, (10) 



J TT 

where in this expression L5 is one of the 0(p A ) couplings of the strong effective chiral 
Lagrangian, and C 4 , C 6 are the Wilson coefficients of the Q4, Qq four-quark operators 
in standard notation. To 0(1/N C ) the scale dependence in Cq cancels with the one 
in <ipil ! >, while C4 is scale-independent (below the charm threshold). The numerical 
result in (|T0"| ) comes from using the chiral limit value f n ~ / = 86 MeV; L5 ~ 1.4 x 10~ 3 
and <Tpip> (lGeV 2 ) = -(0.013 ± 0.003) GeV 3 Q. The Wilson coefficients have been 
evaluated using the perturbative QCD two-loop expressions [0,[T§] restricted to 0(N C ) 
with| Aj^g ~ 300 MeV. The error in (|10"D is partly due to the present error in the 
determination of <i/ji/j>, partly to short-distance uncertainties in C 4j6 . 

The configuration which leads to the b-type coupling is also O(N ). However, 
unlike the case of the c-coupling, we do not know at present how to evaluate b in 
a model independent way, even its leading 0(1/N C ) contribution. There is however 
an important correlation which appears at the order of approximations which we are 
considering: the b-type configuration contributes with opposite signto g 8 and g 27 . This 
correlation of signs is in fact fully respected in the effective action model calculation of 
ref. [fL7[ |, where it is found that b is negative (see also the recent work of refs. 



It is not quite respected in the model of refs. |18|,|19"1|, inspired by the l/A^-expansion; 
but this is due to the fact that in their approach some terms of higher 0(1/N 2 ) have 
also been included. 

A qualitative picture towards the understanding of the underlying physics begins 
to emerge at this simple level of approximations 0(p 2 ) and 0(1/N C ) which we are 
considering. With a= 1 and fixing for example b to the value b ~ —0.52, which is 
the one which follows from the phenomenological determination |g 27 | ~ 0.29, implies 
g g ~ 1.1, still too low compared to the phenomenological number to be explained (which 
once corrected by the enhancement already provided by the 0(p A ) chiral corrections is 
g g ~ 3.6), but in the right direction. 

We are now in the position to bring in as well the discussion of the B^-parameter 
which governs K° — K° mixing. By analogy with the previous analysis of AS = 1 
transitions, the short-distance AS* = 2 Hamiltonian can be visualized as a convolution 
of two s — > d transitions, with two virtual W^-fields being emitted and reabsorbed. 
Each transition results in a flavour matrix factor {Q^ij = ^3^2 times calculable 
short-distance loop functions from the integration of the heavy fields. Therefore, the 
effective 0(p 2 ) AS = 2 Lagrangian has also the structure given in eq. @, but with 
the matrices replaced by Q%2- In this case the configuration c is identically zero 
because (Q32) 2 = while a and b generate the same structure (^ /J ) 23 (£ /i )23. Chiral 
symmetry guarantees that the coefficients a and b appearing in the AS* = 1 and AS = 2 
effective Lagrangians at 0(p 2 ) are the same PDI (once the known short- distance factors 



§ More precisely, we have taken A^T* 00 ~ 400 MeV, which for N c = Nf = 3 corresponds to 
Aj^-SOOMeV. 
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have been appropriately reabsorbed in the global normalizations of £^ s=1 and £^ 5=2 ). 
The expression which emerges for B/^ to lowest 0(p 2 ) in the chiral expansion , but 
taking into account the chiral corrections which bring the chiral limit /-coupling to the 
physical fx, is then 



B 



K 



- (a + b) 



where, at the level of approximation in the l/iV c -expansion that we are considering 
[i.e., 0(1) and 0(1/ N c )], a = 1. Using b = -0.52 as before, one gets B K = 0.36, a 
number which is compatible with the results of the effective action model calculations of 



ref. flT[]] as well as with various phenomeno logical QCD sum rule determinations |2lH23 



Within errors, it is also compatible with the results of the l/A^-approach calculations 
of ref. |p!9j ) but not with the most recent numerical estimates of B/^ obtained by the 



lattice QCD simulations p4]-E£ 



4. In order to get some insight into the underlying QCD dynamics we shall next 
examine the short-distance behaviour of the two-point function correlators 
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in perturbation theory and within the l/iV c -expansion. Here 7i^ s_1 ' 2 (x) denote the 
standard AS* = 1 or 2 four-quark effective Hamiltonians. The spectral functions as- 
sociated to these correlators describe in an inclusive way transitions from the vacuum 
to physical states with total strangeness S = 1 or 2. They have been calculated in 
perturbation theory to next-to-leading logarithmic order in refs. [|TT|,P9"|. The results 
of these calculations give gluonic corrections of rather normal size for the (27l, 1^) 
correlators (i.e., for AS = 2 transitions and AS" = 1 transitions with A J = 3/2) and 
a big enhancement in the (8^, 1r) correlator. The enhancement disappears completely 
when only the large-iV c limit component of the gluonic corrections is retained. 

To simplify the discussion to the essential point let us restrict ourselves to the 
non-penguin operators Q± = Q2 ± Qi and consider the spectral functions associated 
with the C±(fi 2 )Q± terms in the AS = 1 Hamiltonian in the absence of penguin-like 
contributions. The corresponding results from ref. can then be written as follows: 



(12) 
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A very revealing pattern emerges when the coefficients K± of the 0(a s ) corrections 
are expanded in powers of 1/N C as shown above. In the large N c limit K- = K, + 
and the two spectral functions coincide. The 0(1/N C ) corrections to these coefficients 
are enormous, and modify the spectral functions by the same amount but in opposite 

directions: ~Im\I/ gets a large enhancement while -Im\l/ ++ is strongly suppressed. 

Although the higher order 1 / iV c -corrections are smaller than those to next-to-leading 
order, they still have an important overall numerical effect when compared to the exact 
results. This is because in /C+, the alternating signs of the first five terms of the series 
in powers of 1/N C produce a compensating effect, while in K- all the terms have the 
same positive sign which results in an important further enhancement. 

We propose to compare the relative 1 /^-dependence of the spectral functions 
-Im\I/±± calculated in perturbation theory with those obtained to lowest order in xPT 
in the chiral limit and in the l/iV c -expansion. We denote by -Im\I/8,27 the spectral 
functions associated to the effective chiral Lagrangians £g,27 in eqs. @ and @. Then, 
with g 8 = gg + gg , the "equivalent" spectral functions are: 
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We find that the terms of relative 0(1) and 0(1/N C ) in both types of spectral functions 
i.e., those obtained from the effective chiral Lagrangian and those obtained in pertur- 
bation theory, have exactly the same correlation of signs as the one implied by eqs. (H) 
in the limit where a = 1 and in the absence of penguins. 

We shall use this comparison of relative 1 /^-dependence of spectral functions as 
a way to suggest a plausible pattern of the 0(1 /N c ) and 0(1/N^) contributions to the 
couplings a and b. Setting 
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results then in the following equivalence relations: 
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We insist on the fact that these relations are not equalities. The derivation of quanti- 
tative relations would require the use of dispersion relations and a precise knowledge 
of the hadronic spectral functions at intermediate energies which unfortunately is not 
available. The relations above only show the type of 1/N C - corrections in the effective 
couplings which emerge if one assumes that the l/iV c -beliaviour of the short-distance 
correlators is a universal feature of the full hadronic spectral function. The pattern 
suggested by this comparison is nevertheless rather interesting. As expected, the term 
corresponding to (3 is large and negative. It also shows that the 0(1/ N^) term which 
corresponds to a contributes with positive corrections, which tend to cancel in the 
combination a + b. The 0(1/ N^) corrections to b corresponding to f3' have a much 
smaller size. 



5. There are some conclusions we can draw from the previous analyses: 

The phenomenological result |g 27 | ~ 0.29 can be easily digested to lowest order in 
the chiral expansion and to next-to-leading order in the l/iV c -expansion. It requires a 
negative value for the coupling constant b to O(l/N ). This also helps to explain part of 
the AI = 1/2 enhancement, but a quantitative understanding of the phenomenological 
result |g 8 | ~ 3.6 obtained with inclusion of the 0(p A ) chiral corrections is still lacking at 
this level. As already mentioned, a negative b-coupling is a common result of various 
model calculations. However, in order to explain both |g 27 | ~ 0.29 and |g 8 | ~ 3.6 one 
still needs sizable higher 0(1/N^) positive contributions to the a-coupling constant 
which partly compensate in the sum a + b, the large and negative b-coupling which 
is needed to get the AI = 1/2 enhancement. No model so far has been produced 
which shows this convincingly; but it is interesting that both requirements appear to 
be compatible with the pattern of short-distance inclusive calculations discussed above. 

The early chiral symmetry prediction [^(J B/^ ~ 0.35 appears then as a natural 
result within this scenario, but the discrepancy of this prediction with the numerical 
estimates of B^- obtained by the lattice QCD simulations pose a serious puzzle which 
requires further comments on our part. 

If one interprets the large lattice results as chiefly due to the fact that b is a 
very small negative quantity or even a positive one then, from the analysis above, 
it follows that the bulk of the AI = 1/2 enhancement has to come from penguin- 
like configurations i.e., a large and positive value for the c coupling constant^. If 
that is the case we have then to understand why here the l/iV c -expansion, at its first 
non-trivial level, breaks down so dramatically. The QCD perturbative calculation 
of the asymptotic spectral function associated to the penguin Q 6 -operator made in 
ref. [jT7[] shows in fact little difference between the leading result and the one including 
subleading terms in the l/iV c -expansion. There is also another well known problem 
in this case, which is that the predicted value for AI = 3/2 transitions comes out 

' Notice that such large corrections would also imply a large e'/e value. 
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too large: b > 0, results in g 27 > 0.6 i.e., at least a factor of two bigger than the 
phenomenological determination in eq. (||). 

If on the other hand we assume that the estimate b ~ —0.5 is in the right ballpark, 
then the large lattice results for the B^-factor imply that the chiral corrections to 
AS = 2 transitions have to be as large as (9(100%)! Where could such an enormous 
correction come from? The chiral loop corrections to the K° — K° transition amplitude 
have been evaluated by several groups fl9], |3lH33H and it is now known that, once 
the terms which renormalize the / coupling in the chiral limit to the physical fx 
are factorized, the rest of the corrections do not have large chiral logarithmic terms. 
Possible large corrections can then only come from the local 0{p A ) terms of the AS = 2 
effective chiral Lagrangian. The model calculations of these couplings which so far 
have been made |32|,[!3[] give results which are still controversial. These calculations are 
impressive but difficult to interpret. For example, the results of the B^-factor obtained 



in ref. |3B| turn out to be too dependent on the choice of the cut-off which in their 
approach is supposed to separate long- and short-distances contributions. Further 
progress in this direction is indeed possible and hopefully will be made; but we are not 
there yet. 

In the mean time, it seems fair to conclude that there is still, unfortunately, a large 
theoretical uncertainty in our knowledge of the B^- parameter. We do not understand 
the physics behind sufficiently well as yet to restrict the error bars to those of our 
favourite calculation as it is done in many phenomenological analyses of the unitarity 
triangle constraints. It is important to keep in mind that the ultimate purpose of these 
analyses is to test the Standard Model and not some particular QCD estimate of a 
hadronic matrix element. 
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Figure Captions 

Fig. 1: Diagrammatic representation of the three effective couplings in eq. (|6]). The 
solid lines represent quark fields propagating in a gluon background simulated by the 
dotted lines. The (3L^~°P era tors represent the emission and absorption of the virtual 
Vy-field. The restriction to lowest order in the chiral expansion implies that at most 
two £ M insertions are allowed. 

Fig. 2 : Two gluons exchanged from one fermion loop to the other are at least required 
to modify their factorization and this leads to a relative correction of 0(1/N^). 
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